SCALAR PARABOLIC PDE'S AND 

BRAIDS 



Abstract. The comparison principle for scalar second order parabolic 
PDEs on functions u(t, x) admits a topological interpretation: pairs of 
solutions, •) and u 2 (t,-), evolve so as to not increase the inter- 
section number of their graphs. We generalize to the case of multiple 
solutions {u a (t, -)}2=i- By lifting the graphs to Legendrian braids, we 
give a global version of the comparison principle: the curves u a (t, •) 
evolve so as to (weakly) decrease the algebraic length of the braid. 

We define a Morse-type theory on Legendrian braids which we 
demonstrate is useful for detecting stationary and periodic solutions 
to scalar parabolic PDEs. This is done via discretization to a finite 
dimensional system and a suitable Conley index for discrete braids. 

The result is a toolbox of purely topological methods for finding 
invariant sets of scalar parabolic PDEs. We give several examples of 
spatially inhomogeneous systems possessing infinite collections of intri- 
cate stationary and time-periodic solutions. 



We consider the invariant dynamics of one- dimensional second order par- 
abolic equations of the type 



where u is a scalar function of the variables i 6 1 (time) and x G S l = 
M./EZ (periodic boundary conditions in space), and / is a C 1 -function of 
its arguments. 
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1.1. Assumptions. The case of periodic boundary conditions in x 
provides richer dynamics in general than Neumann or Dirichlet bound- 
ary conditions; however, the techniques we introduce are applicable to a 
surprisingly large variety of nonlinear boundary conditions. 

This paper does not deal with the initial value problem, but rather with 
the bounded invariant dynamics: bounded solutions of Eqn. (JTJ) that exist 
for all time t. One distinguishes three types of behaviors which are the 
building blocks of all bounded invariant solutions to Eqn. (JTJ) pjj 114} I22j . 



(i) stationary patterns: u(t,x) = u(x), Vt G M, 

(ii) periodic motions: u(t + T,x) = u(t, x), for some period T > 0, 

(hi) homoclinic /heteroclinic connections: lim^-too u(t, x) = u±(x), where 
u± are stationary or periodic solutions of Eqn. (£[]). 



For the remainder of this paper we impose two natural assumptions on 
Eqn. The first is uniform parabolicity: 

(fl) < A < d w f(x, u, v, w) < A^ 1 , uniformly V (x, u, v, w) G S l x M 3 . 

This condition — that Eqn. (0) grows linearly in u xx — can be relaxed to 
degenerate parabolic equations where the dependence on u xx is as a power 
law, see ^Hl The second hypothesis is a sub- quadratic growth condition on 
the u x term of /: 

(f2) There exist constants C > such that 

\f(x,u,v,w)\ < C(l + H 7 ), uniformly in both x G S 1 and on compact 
intervals in u and w, for some < 7 < 2, 

This will be necessary for regularity and control of derivatives of solution 
curves, cf. This condition is sharp: one can find examples of / with 
quadratic growth in u x for which solutions have singularities in u x . Since 
our topological data are drawn from graphs of u, the bounds on u need to 
imply bounds on u x and u xx . (f2) does just that. 

A third gradient hypothesis will sometimes be assumed: 

(f3) / is exact, i.e., 



f(x, u, u x , u xx ) = a(x, u, u x ) 



- d 
.dx 



d u L - d u L 



(2) 
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for a strictly positive and bounded function a = a(x,u,u x ) and some La- 
grangian L = L(x,u,u x ) satisfying < A < a(x,u, u x ) ■ d\ L(x,u,u x ) < 
A" 1 . 

In this case, we have a gradient system whose stationary solutions are 
critical points of the action J L(x,u,u x )dx over loops of integer period 
in x. This condition holds for a wide variety of systems. In general, 
systems with Neumann or Dirichlet boundary conditions admit a gradient- 
like structure: there exists a Lyapunov function which decreases strictly in 
t along non-stationary orbits. This precludes the existence of nonstationary 
time-periodic solutions. It was shown by Zelenyak [22 that this gradient- 
like hypothesis holds for many nonlinear boundary conditions which are a 
mixture of Dirichlet and Neumann. 



1.2. Lifting the comparison principle. An important property of 
one-dimensional parabolic dynamics is the lap-number principle of Sturm, 
Matano, and Angenent pQ, I16L I20j which, roughly, states that the number 
of nodal regions in x of u(t,x) is a weak Lyapunov function for Eqn. (|T|). 

The lifting of this principle to the simultaneous evolution of pairs of solu- 
tions is extremely fruitful. Consider two solutions v}(t, x) and u 2 (t, x). Any 
tangency between the graphs •) and u 2 (t, •) at time t = t* is removed 
for t = + e (for all small e > 0) so as to strictly decrease the number of 
intersections of the graphs. This holds even for highly degenerate tangen- 
cies of curves [Tj. As shown in the work of Fiedler and Mallet-Paret [10J, 
this comparison principle implies that the dynamics of Eqn. JIJ is weakly 
Morse-Smale (all bounded orbits are either fixed points, periodic orbits, or 
connecting orbits between these), see \14\ I22j . 

The idea behind this paper, following the discrete version of this phenom- 
enon in |12j . is to "lift" the comparison principle from pairs of solutions 
to larger ensembles of solution curves. The local data attached to pairs of 
curves — intersection number — can be lifted to more global data about 
patterns of intersections via the language of topological braid theory. A 
similar theory for geodesies on two dimensional surfaces has been devel- 
oped in [2] , and has served as a guideline for some of the ideas used here. 

Consider a collection u = {u a (t, -)}™ =1 of n > 1 solutions to Eqn. (JTJ), 
where, to obey the periodic boundary conditions in x, {u a (t,0)}™ =1 = 
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{u a (t, 1)}™ =1 as se ts °f points. 1 Instead of thinking of the graphs of 
u a (t, •) as being evolving curves in the (x,u) plane, we take the 1-jet ex- 
tension of each curve and think of it as an evolving curve in (x, u, u x ) 
space. Specifically, for each t, u a (t, ■) : [0, 1] — ► [0, 1] x IR 2 given by 
x i— > (x,u a (x,t),u^(x,t)). As long as these curves do not intersect in their 
3-d representations, we have what topologists call a braid. In particular, 
such a braid is said to be closed (the ends x = and x = 1 are identi- 
fied) and Legendrian (the curves are all tangent to the standard contact 
structure dx^ — x% dx\ = 0). 

As these curves evolve under the PDE, the topological type of the braid 
can change. The topological equivalence class of a closed Legendrian braid 
is the appropriate analogue of the intersection data for pairs of curves. 
Indeed, there is a natural group structure on braids with n strands. We 
argue in a "braid theoretic" version of the comparison principle that the 
algebraic length of a braid given by solutions {u a (t)}™ =1 is a weak Lyapunov 
function for the dynamics of Eqn. (JJJ. 

1.3. Main results. The goal of this paper, following earlier work in 
12 on a discrete version of this problem, is to define an index for closed 
Legendrian braids and to use this as the basis for detecting invariant dy- 
namics of Eqn. ((H). See §21 for definitions and background on the discrete 
version. 

For purposes of detecting invariant dynamics of Eqn. ((TJ, we work with 
braids u relative to some fixed braid v. One thinks of v as a braid for which 
dynamical information is known, namely, that its strands are t-invariant 
solutions to Eqn. (fl)l. the entire set of which respects the periodic boundary 
conditions (individual strands might not: see Fig. [TJ. One thinks of u as 
consisting of "free" strands about which nothing is known with regards to 
dynamical behavior. 

We show that there exists a well-defined homotopy index that maps a 
(closed, Legendrian, relative) braid class represented by {u rel v} to a 
pointed homotopy class of spaces, H(u rel v). This index is at heart a 
Conley index for a suitable configuration space which is isolated thanks to 
the braid-theoretic comparison principle. A coarser homology index sends 
such braids to a polynomial -P T (H) in one variable, r. 



This condition permits solutions with integral period which "wrap" around the 
circle. 
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The main results of this paper are forcing theorems for stationary and 
periodic solutions. 

1.3.1. Stationary solutions. For our main results we restrict to braid 
classes which have two compactness properties: proper and bounded. Roughly 
speaking, a relative braid class {u rel v} is proper if none of the compo- 
nents of u can be collapsed onto u or v. A relative braid class {u rel v} 
is bounded if all strands of u are uniformly bounded with respect to all 
representatives u rel v of the braid class: see §21 

Theorem 1. Let Eqn. satisfy (fl) and (f2) with v a stationary braid. 
If {u rel v} is a bounded proper braid class, then there exists a station- 
ary solution of this braid class if the Euler characteristic ofH, x(H) : = 
P_;l(H), is nonvanishing. If in addition f satisfies (f3) ; then there are at 
least |P r (H)| stationary solutions of this braid class, where \-\ denotes the 
number of nonzero monomials. 

The above theorem is formulated for periodic boundary conditions. In 
the case of other boundary conditions the Zelenyak result implies that 
Eqn. (^) is automatically gradient-like so that the second part of Theorem 
[U is superfluous is those cases. 

Remark 2. With additional knowledge, -P T (H) can reveal more of the 
dynamics. For example, assume for simplicity that the invariant sets are 
known to be hyperbolic and that the strands of u form a single- component 
braid (the graph of u is connected as a subset of S 1 x E). In this set- 
ting, the strong Morse inequalities yield more information on multiplicity 
of solutions. As pointed out before the critical elements of Eqn. (0) are 
equilibrium solutions and periodic orbits. Therefore the Morse relations 
are given by 

« j 

where Q T is a polynomial with non-negative coefficients. The coefficients 
cii count the number of equilibrium solutions of Morse index i, while the 
bj count the number of periodic orbits of Morse index j. If one assumes 
nondegeneracy, then the Morse relations can be used to compute P T (H). 

Remark 3. In the exact case the lower bound on the number of critical 
points can refined even further. For parabolic recurrence relations the 

spectrum of a critical point satisfies A < Ai < A 2 < A 3 < A 4 < A 5 

This ordering has special bearing on non-degenerate critical points with 
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odd index. To be more precise, for a 'topological' non-degenerate critical 
point u with P T (u) = Ar 2k+1 it holds that A = 1. More details of this 
can be found in S J7.3I A direct consequence there are at least as many 
critical points as the sum of the odd Betti numbers of H. If we write 
-P-r(H) = P° dd (H) + P^ ven (H), then our lower bound on the number of 
critical points becomes 

P 1 odd (H) + |P T cvcn (H)|, 
which lies in between |P T (H)| and Pi(H). 

The proof of Theorem Q appears in £0 First, however, we introduce the 
relevant portions of braid theory ($2J), followed by a review (§B1EJ) of the 
discrete braid index constructed in |12j . 

This theory applies to a wide array of inhomogeneous equations. In ^3 we 
show, 

Example 4. The equation 

u t = u xx — - sin 2x u x H cosx „ u(u 2 — 1). (3) 

8 cos x + 

possesses stationary solutions in an infinite number of distinct braid classes. 
As a matter of fact we show that one can embed an Bernoulli shift into 
the stationary equation. 

EXAMPLE 5. For e < 1 and any smooth nonconstant h : S 1 — > (0, 1), the 
equation 

e 2 u t = e 2 u xx + h(x)u(l-u 2 ). (4) 

possesses stationary solutions spanning an infinite collection of braid classes. 
This example was studied by Nakashima |17|, I18j . 

These two examples can be generalized greatly. Theorem EH gives ex- 
tremely broad conditions which force an infinite collection of stationary 
solutions. 

1.3.2. Periodic solutions. We also lay the foundation for using the braid 
index to find time-periodic solutions. For simplicity in the analysis, we 
restrict our attention to equations of the form 

u t = u xx + g(x,u,u x ), (5) 

which trivially satisfies Hypothesis (fl). By assuming Hypothesis (f2) 
(without the w variable), we prove an analogue of Theorem ^ for time- 
periodic solutions of Eqn. (JSJ). As we pointed out before, time-periodic 
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solutions can exist by the grace of the boundary conditions. As the result 
of Zelenyak implies, in most cases a weak version of (f3) holds (gradient- 
like) and the only critical elements are stationary solutions. 

Remark 6. A fundamental class of time-periodic orbits are the so-called 
rotating waves. For an equation which is autonomous in x, one makes the 
rotating wave hypothesis that u(t, x) —U{x — ct), where c is the unknown 
wave speed. Stationary solutions for the resulting equation on £/(£) yield 
rotating waves. Modulo the unknown wave speed — a nonlinear eigenvalue 
problem — Theorem^ now applies. In j3] it was proved that time-periodic 
solutions are necessarily rotating waves for an equation autonomous in x. 
However, in the non-autonomous case, the rotating wave assumption is 
highly restrictive. 

We present a very general technique for finding time-periodic solutions 
without the rotating wave hypothesis. 

Theorem 7. Let Eqn. (jSJ satisfy (f2) with v a stationary braid. Let 
{u rel v} be a bounded proper braid class with u a single- component braid 
and -P T (H) 7^ 0. If the braid class is not stationary for Eqn. ( J5| ). then there 
exists a time-periodic solution in this braid class. 

Remark 8. In certain examples one can find braid classes in which a 
given equation cannot have stationary solutions. Since the only possible 
critical elements in that case are periodic orbits it follows that the Poincare 
polynomial has to be of the form P T (H) = (1 + r)p T (H). The polynomial 
p T (H) gives a lower bound on the number of periodic orbits (in the non- 
degenerate case). The single-component hypothesis on u (namely, that 
the graph of u is connected in S 1 x M) is not crucial. For free strands 
forming multi-component braids u, each component of u will be time- 
periodic. Their periods may not be rationally related, however, leading to 
a quasi-periodic solution in time in the multi-component braid class. 

It was shown in |1J that a singularly perturbed van der Pol equation, 

fH = eu xx + u(l - 5 2 u 2 ) + u x u 2 , 

possesses an arbitrarily large number of rotating waves depending one< 1 
for fixed < 5. We generalize their result: 

Example 9. Consider the equation 



v* = Uxx + ug(u) + u x h(x, u, u x ), 



(6) 
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where the non-linearity is assumed to satisfy (f2), i.e. h has sub-linear 
growth in u x at infinity. Moreover, g and h satisfy the following hypotheses: 

(gl) g(0) > 0, and g has at least one positive and one negative root; 
(g2) h > on {uu x ^ 0}. 

Then this equation possesses time-periodic solutions spanning an infinite 
collection of braid classes. 

We provide details in All of the periodic solutions implied are dy- 
namically unstable. In the most general case (those systems with in- 
dependence), the periodic solutions are not rigid rotating waves. 

2. Braids 

The results of this paper require very little of the extensive theory of braids 
developed by topologists [5]. However, since the definitions motivate our 
constructions, we give a brief tour. 

2.1. Topological braids. A topological braid on n strands is an em- 
bedding (3 : ]j" [0, 1] R 3 of a disjoint union of n copies of [0, 1] into R 3 
such that 

(a) the left endpoints /3(Ui 1°}) are {(°> h 0)}" = i5 

(b) the right endpoints /3(Ui{l}) are 0)}it=i> an d 

(c) j3 is transverse to the planes x\ = constant. 

Two braids are said to be of the same topological braid class if they are 
homotopic in the space of braids: one braid deforms to the other without 
any intersections of the strands. A closed topological braid is obtained if 
one quotients out the range of the braid embeddings via the equivalence 
relation (0,x 2 ,x 3 ) ~ (1, x 2 ,x 3 ) and alters the restrictions (a) and (b) of 
the position of the endpoints to be /?(U™{0}) = /3(U™{1}). Thus, a closed 
braid is a collection of disjoint embedded arcs in [0, 1] x R 2 (with periodic 
boundary conditions in the first variable) which are everywhere transverse 
to the planes x\ =constant. 

In this paper, we restrict attention to those braids whose strands are of the 
form (x, u(x),u x (x)) for < x < 1. These are sometimes called Legendrian 
braids as they are tangent to the canonical contact structure dx 2 — x% dx\. 
No knowledge of Legendrian braid theory is assumed for the remainder of 
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this work, but we will use the term freely to denote those braids lifted from 
graphs. 

2.2. Braid diagrams. The specification of a topological braid class 
(closed or otherwise) may be accomplished unambiguously by a labeled 
projection to the {x\, X2)-plane; a braid diagram. Labeling is done as fol- 
lows: perturb the projected curves slightly so that all strand crossings in 
the projection are transversal and disjoint. Then, mark each crossing via 
(+) or (— ) to indicate whether the crossing is "left over right" or "right 
over left" respectively. 

Since a Legendrian braid is of the form (x, u(x),u x (x)), no such marking of 
crossings in the (x, u) projection are necessary: all crossings have positive 
labels. For the remainder of this paper we will consider only such positive 
braid diagrams. We will analyze parabolic PDEs by working on spaces of 
such braid diagrams. Although Legendrian braids are the right types of 
braids to work with as solutions to Eqn. (cf. the smoothing of initial 
data for heat flow), our discretization techniques will require a more robust 
C° theory for braid diagrams. Thus, we work on spaces of braid diagrams 
with topologically transverse strands: 

Definition 10. The space of closed positive braid diagrams on n strands, 
denoted Q n , is the space of all pairs (u, r) where r G S n is a permutation on 
n elements, andu = {u a (x)}™ =1 is an unordered collection of H 1 -functions 
— strands — satisfying the following conditions: 

(a) Periodicity: u a (l) = u T ^(0) for all a. 

(b) Transversality: for any a ^ a' such that u a (x*) = u a (x*) for 
some x* G [0, 1], it holds that u a (x) — u a (x) has an isolated sign 
change at x = x*. 

Because the strands ofu are unordered, we naturally identify all pairs (u, r) 
and (u, f) satisfying f = ara^ 1 for some permutation o G S n . Henceforth 
we suppress the permutations r from the description of a braid, it being 
understood implicitly. 

The path components of Q n comprise the braid classes of closed positive 
braid diagrams. The braid class of a braid diagram u is denoted by {u}. 
Any braid diagram u with C 1 -strands naturally lifts to a Legendrian braid 
by the 1-jet extension of u a to the curve (x, u a (x), u"(x)). If we allow the 
strands to intersect — disregarding condition (b) of Definition El — we 
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obtain a closure of the space Q n , which we denote fl r 
S n := Q n — Q n defines the singular braid diagrams. 



The 'discriminant' 



2.3. Discrete braid diagrams. From topological braids we have passed 
to braid diagrams in order to describe invariant curves for parabolic PDEs. 
There is one last transformation we must impose: a spatial discretization. 

Definition 11. The space of period al discrete braid diagram on n strands, 
denoted T>^, is the space of all pairs (u, r) where r e S n is a permutation 
on n elements, and u = {-u a }™ =1 is an unordered collection of vectors 



u 



a\d 



strands — satisfying the following conditions: 



(a) Periodicity: = Uq for all a. 

(b) Transversality: for any a ^ a' such that u\ 

(*Ci-<y«vi-<i) <o- 



uf for some i, 
(7) 



As in Definition \1(J\ the permutation r is defined up to conjugacy (since 
the strands are unordered) and will henceforth not be explicitly written. 

The path components of T>2 comprise the discrete braid classes of period 
d. The discrete braid class of a discrete braid diagram u is denoted [u]. If 
we disregarding condition (b) of Definition HH we obtain a closure of the 
space T>2, which we denote The 'discriminant' := T>2 — ^ defines 
the singular discrete braid diagrams of period d. 

Figure ^ summarizes the three types of braids introduced in this section. 




Figure 1 . Three types of braids: a Legendrian topological 
braid [left], its braid diagram [center], and a discrete braid 
diagram [right]. 



2.4. Discretization: back and forth. It is straightforward to pass 
from topological to discrete braids and back again. 

Definition 12. Let u G Q n be a topological closed braid diagram. The 
period-<i discretization of u is defined to be 

disc,(u) = {msc d (u a )} a := {u a (i/d)}«. (8) 
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Conversely, given a discrete braid u G T> r d l , we construct a piecewise-linear 
[PL] topological braid diagram, pl(u) := {pl(m q )} ; where Ph(u a ) is the 
C°-strand given by 

Ph(u a )(x) :=u [d . xi + (d-x- [d ■ x\) (u^ -u [d . xi ). (9) 

The following lemma is left as an exercise. 
Lemma 13. Let u G Q n and v G DJ. 

(1) PL sends the discrete braid class [v] to a well-defined topological 
braid class {pl(v)}. 

(2) Ford sufficiently large, {PL(DISC d (u))} = {u}. 

The second part of this lemma accommodates the obvious fact that braid- 
ing data is lost if the discretization is too coarse. This leads to the following 
definition: 

Definition 14. A discretization period d is admissible for u G Q n if 

{PL(DISC d (u))} = {u}. 

In the next section, we will describe a Morse-Conley topological index for 
pairs of braids which relies on algebraic length of the braid as a Morse 
function. Rather than detail the algebraic structures, we use an equivalent 
geometric formulation of length: 

Definition 15. The length of a topological braid u G Q n , denoted t(u) ; 
is defined to be the total number of intersections in the braid diagram. If 
u G T>2 is a discrete braid, then t(u) := t(PL(u)). 

3. Braid invariants 

We give a concise description of the invariant of [12J for relative discrete 
closed braids. 

3.1. Relative braids. The motivation for the homotopy braid index 
is a forcing theory: given a stationary braid v, does it force some other 
braid u to also be stationary with respect to the dynamics? This necessi- 
tates understanding how the strands of u braid relative to those of v. 

Definition 16. Given v e Vt m , define 

U n REL v := {u G tt n : u U v G tt n+m }. 
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The path components of f2 n REL v, comprise the relative braid classes, 
denoted {u REL v}. In this setting, the braid v is called the skeleton. 

This procedure partitions Q n relative to v: not only are tangencies between 
strands of u illegal, so are tangencies with the strands of v. 

The definitions for discrete relative braids are analogous. 

Definition 17. Given v e Vf, define 

V n d REL v := {u G V n d : u U v G V n d +m }. 

The path components of V d REL v, comprise the relative discrete braid 
classes, denoted [u REL v]. In this setting, the braid v is called the skeleton. 

The operations Disc^ and PL have obvious extensions to relative braids by 
acting on both u and v. 

3.2. Bounded and proper relative braids. 

DEFINITION 18. A relative braid class {u REL v} is called proper if it is 
impossible to find an isotopy u(t) REL v such that u(0) = u, u(t) REL v G 
{u REL v}, fort G [0, 1), and u(l)Uv G X n+m is a diagram where an entire 
component of the braid u(l) has collapsed onto itself, another component 
o/u(l), or a component of v. A discrete relative braid class is proper if it 
is the discretization of a proper topological relative braid class. 

The index we define is based on the topology of a relative braid class. It 
is most convenient to define this on compact spaces; hence the following 
definition. 

Definition 19. A braid class (topological or discrete) is bounded if {u REL v} 
is a bounded set in Q n . 




FIGURE 2. [left] a bounded but improper braid class ; [right] 
a proper, but unbounded braid class. Black strands are fixed, 
grey free. 
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For the remainder of the paper, all braids will be assumed proper and 
bounded unless otherwise stated. 

3.3. The Conley index for braids. Consider a discrete relative 
braid class [u rel v] C T>2 which is bounded and proper. We associate to 
this class a Conley-type index for a class of dynamics on spaces of discrete 
braids. This will become an invariant of topological braids via discretiza- 
tion. 

Denote by N the closure of [u rel v] in the space T>^ REL v. We identify 
an "exit set" on the boundary of N consisting of those relative braids whose 
length i can be decreased by a small perturbation. Let w e dN denote 
a singular braid on the boundary of N and let W be a sufficiently small 
neighborhood of w in T>^ rel v. Then W is sliced by rel v into a 
finite number of connected components representing distinct neighboring 
braid classes, each component having a well-defined braid length i G Z + . 
Define the exit set, N~ , of iV to be those singular braids at which i can 
decrease: 

N~ := CL {w E ON : i is locally maximal on INT(iV)} , (10) 
where CL denotes closure in dN 

Definition 20. The Conley index of a discrete (proper, bounded, relative) 
braid class [u rel v] is defined to be the pointed homotopy class of spaces 

h([u REL v]) = [N/N~] := (N/N~, [N~]) . (11) 

Example 21. Consider the period-2 braid illustrated in Fig. [S[left] pos- 
sessing exactly one free strand with anchor points u± and «2- The anchor 
point in the middle, u\, is free to move vertically between the fixed points 
on the skeleton. At the endpoints, one has a singular braid in E which is 
on the exit set since a slight perturbation sends this singular braid to a 
different braid class with fewer crossings. The end anchor point, u<i, can 
move vertically between the two fixed points on the skeleton. The singular 
boundaries are in this case not on the exit set since pushing w 2 across the 
skeleton increases the number of crossings. 

Since the points u\ and U2 can be moved independently, the configuration 
space N in this case is the product of two compact intervals. The exit set 
N~ consists of those points on dN for which u\ is a boundary point. Thus, 
the homotopy index of this relative braid is [N/N~] ~ S 1 . 
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Figure 3. A period two braid [left], the associated config- 
uration space [center] , and a period six generalization [right] . 



By taking a chain of copies of this skeleton (i.e., taking a cover of the 
spatial domain), one can construct examples with one free strand weaving 
in and out of the fixed strands in such a way as to produce an index with 
homotopy type S k for any k > 0. 

The extension of the Conley index to topological braid diagrams is straight- 
forward: choose an admissible discretization period d, take the Conley in- 
dex of the period-<i discretization, then show that this is independent of d. 
The key step — independence with respect to d — is, unfortunately not 
true. For d sufficiently small, there may be different discrete braid classes 
which define the same topological braid. The information from any one of 
these coarse components is incomplete. The following theorem, which is 
the main result from j!2j . resolves this obstruction. 

Theorem 22 (see [12], Thm. 19 and Prop. 27). For d sufficiently large, 2 
the Conley index ^([DISC^u REL DISC^v]) is independent of d and thus an 
invariant of the topological braid class {u REL v}. 

Definition 23. Given a topological braid class {u REL v} ; define the 
homotopy index to be 

H(u REL v) := /i([DISC d U REL DISC d v]). (12) 

for d sufficiently large. 

For purposes of this paper, the homotopy index is defined with d sufficiently 
large. This is well-defined, but not optimal for doing computations. To 
that end, one can use the more refined formula of j!2j . which computes 
H for any admissible discretization period d via wedge sums: we will not 
require this complication in this paper. 

2 A sufficient though high lower bound is the number of crossings of u with itself 
and with v. 
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For most applications it suffices to use the homological information of the 
index given by its Poincare polynomial 

oo oo 

P T (H) := dim H k(H)r k = dim H *( N > N ') rk - ( 13 ) 

fc=0 k=0 

This also has the pleasant corollary of making the index computable via 
rigorous homology algorithms. 

4. Dynamics and the braid index 

The homotopy braid invariant is defined as a "Conley index." This index 
has significant dynamical content. 

The most basic version of the Conley index has the following ingredients 
jS]: given a continuous flow on a metric space, a subset N is said to be 
an isolating block if all points on dN leave N under the flow in forwards 
and/or backwards time. The Conley index of iV with respect to the flow is 
then the pointed homotopy class [N/N~], where N~ denotes the exit set, 
or points on dN which leave N under the flow in forwards time. Standard 
facts about the index include (1) invariance of the index under continuous 
changes of the flow and the isolating block; and (2) the forcing result: a 
nonzero index implies that the flow has an invariant set in the interior of 
N. In order to implement Conley index theory in combination with braids 
we define the following class of dynamical systems. 

Definition 24. Given d > 0, a parabolic recurrence relation 1Z on Z/cfZ 
is a collection of C 1 -functions IZi : R 3 — > R ; i G Z/c?Z such that for each i, 
d{lZi > and d^R-i > 0. We say that 1Z is exact if there exists a sequence 
of C 2 - generating functions Si such that 

Ki(ui-i,Ui,u i+ i) = d 2 Si-i(ui-i,Ui) + d 1 S i (u i ,u i+1 ) Mi. (14) 
A parabolic recurrence relation (henceforth PRE) defines a vector field on 

|«)=^«_ 1 ,<,< +1 ), (15) 

with all subscript operations interpreted modulo the permutation r: u% +1 = 

u[ . The flow generated by Eqn. (fTHjl is called a parabolic flow on T>2- 
For more details see |12j . Exact PRR's induce a flow which is the gradient 
flowofW(u):= £.£,(<, < +1 ). 
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A parabolic flow acts on discrete braid diagrams in much the same way 
that Eqn. (JTJ) acts on topological braid diagrams. As we have defined 
it in the Conley index for a discrete braid class [u rel v] uses its 
closure N = cl[u REL v] as an isolating block. Indeed, if [u rel v] is a 
bounded, then is a compact set. If [u REL v] is proper, then vector 
field on T>^ rel v induced by 7Z, is transverse to ON, and N is really an 
isolating block for the parabolic flow. The set defined in the previous 
section then is the exit for N. This particular link lies at the heart of the 
theory and follows from the a discrete version of the comparison principle 
jilt H9| - Details of the construction can be found in j!2j . where it is 
shown that the index h([u REL v]) defined via Eqns. (flT)|) and (jllj) is the 
Conley index of any PRR which fixes v. Fig. 0] illustrates the action of a 
parabolic flow on braids. 




FIGURE 4. A parabolic flow on a (bounded and proper) 
braid class is transverse to the boundary faces, making the 
braid class into an isolating block. The local linking of 
strands decreases strictly along the flow lines at a singular 
braid u. 

In |12j it furthermore is shown that certain Morse inequalities hold for 
stationary solutions of Eqn. (jl5j) . The Morse inequalities also provide 
information about the periodic orbits. This is due to the fact that for 
parabolic systems the set of bounded solutions consists only of stationary 
points, periodic orbits, and connections between them. 

Theorem ^ is an extension of the following results for parabolic lattice 
systems. 
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Theorem 25. [12\ Let 71 be a parabolic recurrence relation. The induced 
flow on a bounded proper discrete braid class [u rel v], where v is a 
stationary skeleton, has an invariant solution within the class [u rel v] if 
the Conley index h = h([u REL v]) is nonzero. Furthermore: 

(1) If the Euler characteristic x(h) ^ then there exist stationary 
solutions of braid class [u rel v] . 

(2) If 7Z is exact, then the number of stationary solutions of braid 
class [u rel v] is bounded below by \P T (h)\, the number of nonzero 
monomials of the Poincare polynomial of the index. 

If a proper bounded braid class [u rel v] contains no stationary braids 
for a particular recurrence relation TZ, then h(u REL v) ^ forces peri- 
odic solutions of Eqn. (|15|). i.e. the components of u are periodic. If the 
system is non-degenerate the number of orbits is given by P\(h)/2. As a 
consequence in this case P T (h) is divisible by 1 + r and 71 is not exact. 
Note that for d large enough the topological information is contained in 
the invariant H for the topological braid class {u rel v}. 



5. Examples: stationary solutions 

The following examples all satisfy Hypotheses (fl) and (f2). 

Example 26. Consider the following family of spatially inhomogeneous 
Allen-Cahn equations studied by Nakashima |17|, I18j : 

e 2 u t = e 2 u xx + h(x)u(l - u 2 ), (16) 

where h : S 1 — > (0, 1) is not a constant. Clearly this equation has stationary 
solutions u — 0, ±1 and is exact with Lagrangian 

L= l -e 2 u 2 x -\(h{x)u 2 {2-u 2 )). 

According to |17j . for any N > 0, there exists an > so that for all 
< e < e^r, there exist at least two stationary solutions which intersect 
u = exactly iV times. (The cited works impose Neumann boundary 
conditions: it is a simple generalization to periodic boundary conditions.) 

Via Theorem |2H1 we have that for any such h and any small e, this equation 
admits an infinite collection of stationary periodic curves; furthermore, 
there is a lower bound of iV on the number of 1-periodic solutions. 
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Example 27. Consider the following equation 



u t = u xx — ^ sin 2x u x H """"" 3 u(u 2 - 1), (17) 

o cos x + -7= 



5 . _ cosx 
with x eS 1 = M/2ttZ. 

Eqn. © is a weighted exact system with Lagrangian 



, ,«•>,. I ' » cosx ill 2 - 1) 

L = e ie cosZx -ui 



2 x cos x + -4= 4 
where by "weighted exact" we mean (cf. Eqn. ((2))) 

— ros2r I d 8L 8L 



(19) 



One checks easily that there are stationary solutions u = ±1 and u± = 
±| (-\/5cosx+ l), as in Fig. These curves comprise a skeleton v = 
{ — 1, m_, m + , +1} which can be discretized to yield the skeleton of Exam- 
ple |^ From the computation of the index there, this skeleton forces a 
stationary solution of the braid class indicated in Fig. EJleft]: of course, 
this is detecting the obvious stationary solution u = 0. 

What is more interesting is the fact that one can take periodic extensions 
of the skeleton and add free strands in a manner which makes the rela- 
tive braid spatially non-periodic. Let us describe a family of proper and 
bounded relative braid classes. Let v n be the n-fold periodic extension of 
v on [0,n — 1] and consider a single free strand — 1 < u(x) < 1 that links 
with v n as follows: on each interval [k, k + 1], k = 0...n — 1, we choose one 
of three possibilities: 



(a) l{u,uJ) = and i(u,u + ) = 2, 

(b) i(u,U-) = 2 and i(u,u + ) = 2, or 

(c) i(u,U-) = 2 and i(u,u + ) = 0. 



Define a symbol sequence a = (ai...a n ), where Oi G {a, b, c}. Every symbol 
sequence except for a = (a... a) and a = (c.c), defines a proper and relative 
braid class {u^ rel v}. 
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To compute the invariant, we discretize. Choose the discretization d = 2n 
on [0,n]. 3 Fig. EJright] shows an example. In JJJthe index was computed: 

H(u CT REL v) = /l(DISC 2n U (J REL DISC 2n V n ) ~ S k , 

where k = G a}. Therefore, P r (H) = r k . 

The Morse inequalities now imply that for each n > there exist at least 
3™ — 2 different stationary solutions. This information can be used again to 
prove that the time-27r map of the stationary equation has positive entropy. 




Figure 5. The skeleton of stationary solutions for Eqn. (J3J) 
forces an infinite collection of additional solutions which 
grows exponentially in the number of strands employed. 

Example 28. The following class of examples is very general and includes 
Example I2IH as a special case. One says that Eqn. (JTJ) is dissipative if 

uf(x, u, 0, 0) — ► — oo as \u\ — ► +oo (20) 

uniformly in x G S 1 . 

Theorem 29. Let f be dissipative and satisfy (fl)-(f2). J/v is a nontriv- 
ially braided stationary skeleton (i.e., t(v) ^ 0), then there are infinitely 
many braid classes represented as stationary solutions to Eqn. (jT|). More- 
over, the number of braid classes for which u consists of just one strand is 
bounded from below by 2\l/2], where i is the maximal number of intersec- 
tions between two strands of v. 

Proof. Given the assumptions one can find c + > and c_ < such 
that ±f(x, c±, 0, 0) > 0, and 

c_ < v a (x) < c + , 

3 This is admissible for the skeleton v" and is large enough to yield the correct index 
computation. 
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for all strands v a in v. Using discrete enclosure via sub/super solutions, 
Lemma in Appendix C yields solutions u + and w_ such that 

c_ < u-(x) < v a (x) < u+(x) < c+, 

for all a. Assume without loss of generality that all strands in v are 
1-periodic (if not, one can take an appropriate covering of v). For the 
sake of convenience we may assume that x G S 1 = K./Z. Select two in- 
tersecting strands which form the braid w = {v ai ,v a2 }, and set t(w) = 
^{intersections between v ai andf Q2 }. Consider the skeleton z = {y,-,v ai ,V 
and a free strand u(x) - - with u(x + 1) = u(x) - - that links with z 
as follows: (1) u~(x) < u(x) < u+(x), (2) for some k > 0, i(u,v ai ) = 
i(u,v a2 ) = 2k < l(w). These two hypotheses describe the relative braid 
class {u rel z}, which clearly is a proper and bounded class and therefore 
has a well-defined homotopy braid index H. The index H is an invariant of 
the braid class and it can be computed for instance by studying a specific 
system of which all solutions are known. 

Consider the equation e 2 u xx +u — u 3 = 0. If we choose e = (7r(t(w) + 1))~ , 
then there exists a periodic solution v\{x) with period T = 2/t(w). Define 
V2{x) = V\ (x— (l/t(w))j ; then if we consider v\ and v% on the interval [0, 1], 
it follows that l(v%,V2) = t(w). The skeleton z' = {— 1, V\, v 2 , +1} is now 
topologically equivalent to w. Moreover, the equation e 2 u xx +u— u 3 = has 
a unique solution u which has the right linking properties with the skeleton 
z': —1 < u{x) < 1, and i(u, v%) = l(u, v 2 ) = L{u,v ai ) = i(u,v a2 ) < t(w). 
Therefore, {u rel z'} and {u rel z} are topologically equivalent. As in 
|12[ |2j the invariant set Inv({u rel z'}) of the equation 

u t = e 2 u xx + u-u 3 , 

is given by Inv({u rel z'}) = {u(x + <p) | <p & R}, which represents a 
hyperbolic circle of stationary strands. Its unstable manifold has dimension 
i(u, Vx) = 2k and therefore its Morse polynomial is given by r 2fe_1 (l + r). 
Since this captures the entire invariant, it follows that P T (H) = r 2fc_1 (l + 
r): see also jH] for details. 

From the invariant H and Theorem ^ we deduce that if Eqn. (0) is dissipa- 
tive and exact it has at least pairs of 1-periodic solutions. One finds 

infinitely many stationary braids by allowing periods 2n. Indeed, take the 
periodic extension w 2 ™. Then for any k satisfying 2k < i(w 2n ) = 2rw(w) 
we find a 2n-periodic solution. By projecting this to the interval [0, 1] we 
obtain a multi-strand stationary braid for Eqn. ((TJ. As a matter of fact 
for each pair p, q, with q < p and gcd(p, q) = 1, there exists at least two 
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distinct periodic solutions v} p q and u 2 by setting k = i(w)q and n = p. 
This infinity of solutions enshrouds the set Q fl (0, 1). □ 

6. Examples: time periodic solutions 

This is a longer example of a very general forcing result for time-periodic 
solutions. 

Example 30. Consider equations of the following type 

Ut = u xx + ug{u) + u x h(x,u,u x ), x G R/Z, (21) 

where the non-linearity is assumed to satisfy (f2), i.e. h has sub-linear 
growth in u x at infinity. Moreover, assume that g and h satisfy the hy- 
potheses: 

(gl) g(0) > 0, and g has at least one positive and one negative root; 
(g2) /i > on {uu x ^ 0}. 

Theorem 31. Under the hypotheses above Eqn. (|2Tj) possesses an infinite 
collection of time-periodic solutions all with different braid classes. 

Proof. Consider first the perturbed equation, 

u t = u xx + ug(u) + a e u x h(x,u,u x ), (22) 

where a e = for ^u 2 + u 2 x G [0, e] and a e = 1 for ^Ju 2 + u 2 x > 2e. For 
e > Eqn. (}2*2"|) has small stationary solutions u e (x) which oscillate about 
u = 0. We can choose this u e and a relatively prime pair of integers p, q G N 
such that u e (x + p) = u e (x) and g{fS)/2n < q/p is arbitrarily close to 
q/p. The integer q represents the number of times the oscillation fits on 
the interval [0,p\. 

We use (gl) to build a skeleton for Eqn. (j2^|l . Let a + and a_ denote pos- 
itive and negative roots of g, and consider the skeleton v = {v 1 , v 2 , v 3 , v 4 } 
on M/pZ with v 1 (x) = a_, v 2 (x) = a + , v 3 (x) = u e (x), and v 4 (x) = u t (x — 
p/2q). Clearly l(v 3 , v A ) = 2q. Define the relative braid class {u rel v} 
as follows; u = {u} is a (1-strand) braid satisfying a_ < u(x) < a + and 
t(u, v 3 ) = l{u, v 4 ) = 2r < 2q. This braid class is proper and bounded, and 
its homotopy invariant H was computed in the previous section: 

P r (H) = r 2 ^ 1 (l + r). 
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We claim that for < e < 1 there are no stationary solutions in {u REL v}. 
Suppose that u is stationary. One checks that the function 

H{u,u x ) := -u 2 x + u I g(s)ds — / / g(r)drds 
2 Jo Jo Jo 



has derivative 



~rH = —a e u x h(x,u,u x ). 
dx 



This term is nonpositive by (g2) and not identically zero by the fact that 
u cannot be close to a constant (thanks to the intersection numbers). The 
periodic boundary condition leads to the desired contradiction. 

Since u is a 1-strand braid it follows from Theorem [7| that {u rel v} 
contains a t-periodic solution to Eqn. ((22J). By lifting the equation to 
the interval [0,fcp], k e N, we obtain different periodic solutions for each 
r < kq, which shows that there are t-periodic solutions for infinitely many 
different braid classes: see [12, Lem. 43] for details. What remains is to 
show that these periodic solutions to Eqn. (|2*2*j) persist in the limit e — > 0. 
We need to show that the limits obtained are not equal to the zero solution. 
We use an argument similar to that of Angenent [5] . 



Linearize Eqn. (|22jl around u = 0. This leads to the linear operator L = 
d 2 /dx 2 +g(0) on L 2 (M./pZ). The spectrum of L is given by the eigenvalues 
\ n = -4ir 2 n 2 /p 2 + 0(0), for n = 0, 1, .... Since y/g(0)/2n < q/p it holds 
that A„ > for all n < q, and A n < for n > q. This yields a (spectral 
decomposition) splitting of L = L + + L_. The evolution on the set I = 
{ip | t,(ip, 0) = 2r < 2q} is then dominated by the linear operator L for 
HV'IU 2 small. Therefore, the function B(ip) = ~ (ip, L + ip) L 2 satisfies 

j t B{iP) = {L+iP,L + iP) L 2 + o{U\\ L 2). (23) 

for all ip G I. For u e (t,x) a periodic solution with sufficiently small L 2 
norm, Eqn. (}2"3*|) implies that -j^B > 0, a contradiction of periodicity. Thus 
we conclude that the e — > limits do not collapse to zero. □ 

REMARK 32. The form of Eqn. (J2"T|) is not the most general form possible. 
Certainly, having h strictly negative on {uu x ^ 0} is also permissible. With 
work, the diligent reader may verify that allowing the u xx term to vary as 
per (fl) does not change the nature of the results. 



SCALAR PARABOLIC PDE'S AND BRAIDS 23 

7. Proofs: Forcing stationary solutions 

7.1. Discretization of the equation. From hypothesis (fl) we ob- 
tain an estimate for / of the form 

f(x, u, v, 0) + a_{w)w < f(x, u, v, w) < f(x, u, v, 0) + a + (w)w, (24) 

for all x G [0, 1], and u,v,w G R, where a_(s) = A -1 for s < 0, a_(s) = A 
for s > 0, and a+(s) = A for s < 0, a + (s) = A -1 for s > 0. 

Consider a braid u of n strands. For the remainder of this section, we work 
with individual strands u = u a , suppressing the superscripts for notational 
aesthetics. 

We discretize Eqn. (0) in the standard manner. Choose a step size 1/d, for 
d G N, and define := u(i/d). We approximate the first derivative u x (i/d) 
by Aui := d{u i+ i — Ui) and the second derivative u xx (i/d) by A 2 Ui : = 
d 2 (w i+ i - 2ui + Ui-i). 

Lemma 33. Let u be a stationary braid for Eqn. (|Tj) . t/ien 

ei(d) := / w is An,, A V)) — > 0, (25) 
as d — > oo uniformly in i. In particular, \€i(d)\ < C/d. 

Proof. From Appendix A it follows that each strand u of a stationary 
solution to Eqn. is C 3 . A Taylor expansion yields 

Aui — u x = d ■ (u((i + 1) / d) — Ui(i / d)) — u x (i / d) 

= d ■ R\ /d (i/d) = ^u xx (y)/d, for some y 

A 2 Ui -u xx = d 2 ■ (u((i + l)/d)-2ui(i/d)+u((i-l)/d)) -u xx (i/d) 

= d 2 ■ R 3 1/d (i/d) = ^u xxx {y)/d. for some y 

For Xq :— i/d it therefore holds that 

\Au d . X0 - u x (x )\ < C/d, lA 2 ^.^ - u xx (x )\ < C/d, 

with C independent of xq. Since / is C l the desired result follows. A more 
detailed asymptotic expansion for is obtained as follows: 

€i(d) = f(i/d,Ui,Aui,A 2 Ui)-f(i/d,u(i/d),u x (i/d),u xx (i/d)) 

= d u J{i/d,u{i/d),u x {i/d),u xx {i/d))[Aui - u x (i/d)] 

+d Ux J(i/d,u(i/d),u x (i/d),u xx (i/d)) [A 2 Ui - u xx (i/d)] + R 2 . 
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From the weak form of Taylor's Theorem the remainder term R2 satisfies 
\R 2 \ = o(l /d). Combining this with the estimates obtained above we derive 
that \ei(d)\ < C/d, thus completing the proof. □ 

The next step is to ensure that the (i-point discretization of u is in fact a 
solution of an appropriate parabolic recurrence relation. 

Lemma 34. Let u e Q n , and let d be an admissible discretization for u. 
Then for any sequence {e"} with i — 0, .., d and a — 1, .., n, there exists a 
parabolic recurrence relation £f satisfying 

^«_ l3 <,< +1 ) = -e« 

where uf = u a (i/d). In addition, \Sf(r, s,t)\ < C maxj iQ , \ef\ for all \r\, \s\, \t\ < 
2maxj )Q \uf\ and some uniform constant C depending only on u. 

Proof. This proof is a straightforward extension of |12[ Lemma 55], 
in which ef = 0. □ 



Let v e Q m be stationary for Eqn. (0) and let {u rel v} be a bounded 
proper braid class with d a sufficiently large discretization period. We now 
construct a parabolic recurrence relation for which the discrete skeleton 
DisCrfV is stationary. Combining the Lemmas ESI and 02J the recurrence 
relation defined by 

lZf(ui-i,Ui,Ui +1 ) :=f(i/d,u u Au i ,A 2 u, i )+£f(u i ^ 1 ,Ui,u i+1 ), (26) 

has DisCdV as a stationary solution. The above construction works for 
any d' > d. For a given d the recurrence relation TZf is considered on the 
compact set cl[u rel v], which implies that \ui\ < 2maXj iQ , \vi\. To verify 
the parabolicity of IZf , we compute the derivatives. From hypothesis (fl) 
and the parabolicity of £f, we obtain: 

d{R* = ■ d 2 + diSf > X ■ d 2 > 0. (27) 

OUr, 



"XX 



Furthermore, 



B ,- Ri= BL.t + *L. i+a ft, 
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which does not yet prove parabolicity since no estimates for d Ux f are given. 
However, utilizing hypothesis (f2), we have that 

TZi (ui- 1} Ui,Ui + i) = f(i/d,Ui,Aui,A 2 Ui)+Si(ui^i,Ui,u i+1 ) 

> f(i/d, u h Au u 0) + a_ (A 2 Ul )A 2 Ui + Sf 

> —C — C\Au i \ y + \A 2 Ui 

> -C -CcT + \d 2 (u i+1 -2u i + u i - 1 ), 

which shows that TZf is an increasing function of u i+ i, provided that d 
is large enough. This relies on the fact that the braid class is bounded. 
Periodicity of / implies that Rf +d = TZf. Initially Hf(r, s,t) is defined for 
\r\, \s\, \t\ < 2maxj jQ , \vi\. It is clear that TZf can easily be extended to a 
parabolic recurrence relation on all of R 3 . 

7.2. Convergence to a stationary solution. Choose d* large enough 
such that TZf is parabolic for all d > d*. Let {uf' d } be a sequence of braids 
which are solutions of 

1Zf(ui-i,Ui,u i+ i) = f(i/d,Ui,Aui,A 2 Ui) + £f(ui-x, v-i, Ui+i) = 0, (28) 

and which satisfy the uniform estimate < C for all d. For notational 
simplicity, we omit the discretization period d and write instead of u*' 
in what follows. The discretization index will be clear from the range of 
the index i. 

Lemma 35. Let {iii} satisfy lZf(ui-i,Ui, Uj+i) = and \m\ < C as d — > oo. 
Then 

d ^ d 

J2d^- C ; X>k+i-^| 2 <C> (29) 

8=0 1=0 

with C independent of d. 

Proof. For each strand a it holds that either Ui + d = Ui, or Ui + kd = Ui 
for some k. Since there are only finitely many strands, the constant k 
can be chosen uniformly for all a. Therefore we assume without loss of 
generality that the first equality holds. The first estimate immediately 
follows from the uniform bound on itj. 

From Eqn. (|28jl it follows that 

f(i/d, Ui, Am, A 2 Ui) = -£f(ui-t, u h u i+1 ). 
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Multiply the above equation by Ui/d, then from Eqn. (J2%|) it follows that 

-A 2 m ■ -Ui < —f(i/d,Ui,Aui,0)-Ui + —Sf{ui^u h u i+1 )-u h 
a a + a a + a 

for Ui > 0, and 

-A 2 iii ■ -Mi < +—f(i/d,Ui,Aui,0)~Ui + — w i; 

da- da- d 

for Ui < 0. From the periodic boundary conditions it follows that 

d d 



d j d 



i=0 i=0 



2 



Combining the above estimates and using (f2) we obtain 

d d ^ ^ 

^d-\u i+1 -Ui\ 2 < ^2^-\f(i/d,Ui,Aui,0)\\ui\ + -\£i\\ui\^ 



-A,„ A 

I 

i=0 i=0 

|2 



< C/d + J2^[Ce + ed 2 \u l+1 - 

i=0 

d 

< C e + eC d ■ \Ui+\ — Ui\ 2 , for any e > 0. 



i=0 

Choosing e small enough yields the second estimate. □ 

Define <J) d := Ph({ Ui }). Then WMh < Etojl^^ and \\£M\h = 
Y^if=o d ' \ u i+i ~ u i\ 2 - Due to the uniform estimates we obtain the Sobolev 
bound < C, with C independent of d. Therefore, (pd„ converges 

to some function u G C°([0, 1]), as d n — > 00. 

Lemma 36. Let {u^ satisfy 1Zf{ui-i,Ui,u i+ i) = and \ui\ < C as d — > 00. 
Then 

d 

J^d- \Aui - Am-^ < C, (30) 
with C independent of d. 

Proof. As in the proof of Lemma l3*K| we have 

A 2 Ui < — —f(i/d,Ui,Aui,0) - —£f(ui-i,Ui,u i+1 ), 
a- a- 
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for A 2 Ui > 0, and 

-A 2 u { < —f(i/d,u h Aui, Q) + — £f (u^, m, u i+ i), 
a + a + 

for A 2 Ui < 0. Combining these estimates with (f2) we obtain 

d- \Am- Auti-tl < C + ClAuil 1 . (31) 

Therefore 

d 



Y^d^~ l \A Ui - A Ul ^\ 2h < C + C^^Aui 

i=0 i=0 



2 



< C by Lemma E3 
which is the desired estimate. □ 

Set ^ d := PL({A«i}). Then ||VJ£ 2 < ^to W^? < C i and WtMWv, = 
Yli=o ' | Aitj + i — Aui\ 2 ^ < C. This implies that HV'dllifi.a/T < C, inde- 
pendent of d. Therefore there exists a subsequence ip dn converging to some 
function v G C°([0, 1]). 

Lemma 37. Let {%} satisfy TZf(ui-±, Ui, u i+ i) = 0, and \ui\ < C as d — ► oo ; 
then 

\A Ui \ < C, |A 2 Ui| < C, 
with C independent of d. 

Proof. The first estimate follows from the fact that || V'rf || c° < C, 
hence \Av,i\ < C. For the second estimate we use Eqn. (|31j) . The uniform 
bound on Aui then yields a uniform bound on A 2 Ui. □ 

Finally, we require an estimate on A 3 Ui = d ■ (A 2 Uj + i — A 2 Ui). 

Lemma 38. Let {ui} satisfy TZf (««_!, Ui, ttj+i) = and \ui\ < C as d — > oo. 
Then 

\A 3 Ui \ = d ■ \A 2 u i+l - A 2 Ui \ < C, (32) 
with C independent of d. 

Proof. Since TZ d d+1 - TZ d = it follows from the definition of TZf that 

f((i + l)/d, u i+ i, Au i+1 , A V+i) - f{i/ d, u h Au h A 2 Ui) 
= -£t +1 {ui,u i+1 ,u i+2 ) +£i(u i _ l ,u i ,u i+ x). 
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Using Taylor's theorem we obtain that 
d x f(i/d,u i ,Au i ,A 2 u i )^ + d u f(i/d,Ui, A%, A 2 Ui)(u i+1 - Ui) 

+ d Ux f(i/d,u i ,Au i ,A 2 u i )(Au i+1 - Aui) 
+ 9 Uxx f(i/ d, Ui, Au^ A 2 Ui) (A 2 u i+ i - A 2 u { ) 
= -{Sf +1 - St) - R 2 {i/d, u h A Ul , A\i). 

For A 3 Ui this implies 

d Ux J{i/d,Ui,Au i ,A 2 u i )A 3 u i = -d x f{i/d,Ui,Aui,A 2 Ui) 

- d u f(i/d,u i ,Au i ,A 2 u l )Au i 

- d Ux f(i/d,Ui,Aui,A 2 Ui)A 2 Ui 

- AS? - R 2 (i/d,Ui,Aui,A 2 Ui)d. 

By Lemma IH7I the right hand side is uniformly bounded in d. Using (fl) 
then yields the desired estimate on A 3 Wj. □ 

Define Xd '■= PL({A 2 Wj}). From LemmalHHlwe then derive that H^XrflU 00 — 
C. Therefore Xd„ converges to some limit function w in C°([0, 1]). 

From Lemmas HUE ESI and 1381 it follows that the functions <fid, ipd and Xd 
converge to function u, v and w respectively, with the anchor points being 
solutions of TZf = 0. 

The following lemma relates discretized braids to stationary braids in 
{u rel v}. 

Lemma 39. Let u,v,w e C°([0, 1]) and let {uf}f =0 be sequences whose PL 
interpolations satisfy 

PL(uf) -> u, PL(Awf) -> v, PL(A 2 uf) -> w, 

in C°([0,1]) as d -> oo. If Tlfiuf^uf,^^) = and |A 3 «j| < C, toen 
u G C 2 ([0, 1]) ; m x = f, andu xx = w satisfying f(x,u,u x ,u xx ) = pointwise 
on [0,1]. 

Proof. We start with the estimate \<fi' d — ipd\ < 2 max o<;<d I A 2 it;| — > 
uniformly as d — > 00. This implies that ipd v in C°([0, 1]). The 
same estimate holds for \ip' d — Xd\ < h max <i<d | A 3 Ui\ — > uniformly 
as d — > 00. Hence we deduce that Xrf ~~ ► w m C°([0> 1])- From the 
definition of derivatives it now follows that HZV^u — v —* 0, and 
ll-Di/df — HU°° — > 0; thus v = u x and w = u xx . From Lemma I3H1 we 
deduce that f(x, u, u x , u xx ) =0. □ 
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Note that Lemma PHI in Appendix A implies further that u G C 3 ([0, 1]). 

7.3. Proof of Theorem QJ Given P_iH ^ 0, the existence of a single 
stationary solution is argued as follows. Choose large enough. Then 
from Theorem ESI it follows that Eqn. ()28)1 has a discrete braid solution 
{u" ,d } for all d> d*. The boundedness of the braid class implies that the 
sequence } satisfies 

\v% ,d \ < C, V i, a, and V d > d*. 

Lemmas imply that as d — > oo one obtains a stationary braid u = 

{u a } whose strands u a satisfy the equation f(x, u a , u°, u" x ) = 0. Since the 
skeletal strands DlsCd(f /3 ) converge to v@ by construction and the pairwise 
intersection numbers are the same for all d, we have in the limit a solution 
to Eqn. (TjQ) in the correct braid class. 

It remains to determine multiplicity in the case of (f3). The difficulty lies 
in dealing with degenerate critical points: one proceeds using the standard 
tools of critical groups and Gromoll-Meyer pairs. We refer the interested 
reader to |7] for detailed definitions. For the remainder of the proof, we will 
characterize Morse data of critical points u via the Poincare polynomial 
-P-r(u). For a nondegenerate critical point, this is a polynomial of the form 
P T (u) = r^ u \ where fi is the Morse index. For degenerate critical points, 
P T is defined via certain homology groups pp. 

In the gradient case one has the action A on the space Q n defined as follows: 

A(u) = J2 L(x,u a ,u^)dx, uGlF, 
«=i Jo 

and the discretized action on defined by 

n d n d 

A d (u) = £E L G> <, A<) + E E a <«' A <)' u e ^ 

q=1 i=0 a=l i=0 

where the a, are small perturbations guaranteeing that Disc^v is a critical 
skeleton for each d > (i*. 4 These can be constructed as in Lemmas l33ll34l 
so as to satisfy the same estimates. It follows immediately from Eqn. (J2J) 
that Kf = -d Ui A d . 

Assume without loss of generality that A has finitely many critical points 
Uj REL v so that all critical points are isolated. We have shown earlier in 



We omit the superscript d in the notation for uf. 
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this section that as d —>■ oo, critical points of A4 converge to a critical point 
of A. We will factor this convergence through a sequence of nondegenerate 
Morse functionals in order to extract forcing data. 

One may perturb A on a neighborhood of the critical points to A £ which is 
Morse on the braid class {u rel v}. Next, discretize A e to yield functionals 
A d . Our convergence results imply that A e d is Morse for d sufficiently large. 
Indeed, if {u d } is a sequence of critical points of A € d , then PL(u d ) converges 
in Q n to a critical point of A. The same holds for the eigenfunctions and 
eigenvalues of the linearized functional, which implies that A e d is Morse for 
d large enough. Uniform estimates on the remainder terms of A e and A d 
then yield uniformity in the distance between the critical points of A d for 
all d large. To be more precise dis%^(u^, u^,) > 5 t > for any pair of 
critical points u^, u^, of A d . 

Let Bi be small isolating neighborhood of the critical points u« of A. For 
e > sufficiently small all critical points of A e are contained in the neigh- 
borhoods Bi. We can group together the critical points of A% in asso- 
ciated neighborhoods Bf in T> d . Since the critical points form a Morse 
decomposition for A we obtain the following Morse inequalities from j!2l 



where Q T has nonnegative coefficients. Due to the uniform separation of 
the critical points as d — > 00 Eqn. (|33|) also holds in the limit for the 
functional A\ i.e. P T (A) = P T (H) + (1 + r)Q T . 

Lemma EH1 below shows that for a given braid class, each critical point u, 
of A has Poincare polynomial of the form P T (uj) = AiT Pi . By Lemma l4ll 
following, we can find Morse approximations A e whose Poincare polynomial 
is exactly the same, i.e. P T (Bi) = AiT Vi . Substituting the latter into then 
Morse inequalities for A e the proves that the number of neighborhoods 
Bi is bounded from below by the number of monomials in P T (H) — i.e. 



Lemma 40. Given u an isolated critical point of A, the Poincare polyno- 
mial is of the form P T (u) = Ar p for some AeN and p > 0. 



§7] 




(33) 



Pr(H)|. 



□ 



Proof. In the case of a braid class with a single free strand, the conclu- 
sion follows from a result of Dancer |9] : since A is a first order Lagrangian 
of a scalar variable, a degenerate critical strand has nullity at most two. 
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In the case of braids with multiple free strands, the proof becomes some- 
what more delicate. By considering the appropriate covering we obtain 
an uncoupled system of equations for the components of the braid u. The 
critical groups of the braid class are precisely the tensor product of the crit- 
ical groups of the individual components (see Theorem 5.5 of |7J). Thus, 
the Poincare polynomials multiply, and the result follows from the single- 
strand case. □ 

Lemma 41. Given A having finitely many critical points Uj e Bi with 
P T (uj) = AiT Pi , there exists a C 2 -small perturbation of A with support in 
Uj-Bj to a Morse functional A e having exactly A4 critical points in B i; each 
with Morse index pi. 

Proof. We consider each degenerate critical point separately. For 
each degenerate critical point, the data in its critical groups comes from 
a 2-dimensional 'center' set W given by the Gromoll-Meyer version of the 
Morse Lemma [13J: all the non-hyperbolicity of dA is manifested on W. 

Consider ^4|vy : M 2 — *■ M with coordinates chosen so that there is a degen- 
erate critical point at the origin having P T = A4T. The statement of the 
lemma now becomes the claim that there exists a perturbation of A\w to 
a function on R 2 which has A4 critical points of Morse index one. This 
follows from choosing a small disc D at the origin which is an isolating 
neighborhood for VA. (This is possible via a result of |21j .) This implies 
that VA is transverse in/out of 3D on an alternating sequence of 2Ai + 2 
arcs as in Fig. deleft] . 

One may then set up analytic coordinates on D and write out an explicit 
Morse function with Aj saddle points. A less explicit method is to note 
that a linear chain of Aj saddles — as in Fig. fright] - possesses an 
isolating neighborhood whose boundary is combinatorially equivalent to 
that of the disc D: for D small, mapping this chain of saddles to D yields 
the appropriate perturbation of A. □ 



8. Proofs: Forcing periodic solutions 

In this section, we provide details of the forcing arguments in the case of 
non-stationary solutions. The technique is philosophically the same as for 
stationary solutions: discretize, apply the Morse-theoretic results of [12J, 
then prove convergence to solutions of Eqn. ((TJ. However, the requisite 
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Figure 6. [left] An isolating neighborhood of the critical 
point of A\w', [right] Replace it with a chain of nondegenerate 
saddles. 



estimates are more involved in the time-periodic case. For this reason, we 
present the proofs for the normalized equation, 

u t = u xx + g(x,u,u x ), (34) 

noting that the general case of Eqn. ([TJ is valid, though messier. 

Appendix|B]details a regularity result for non-stationary solutions to Eqn. (p£ij). 



8.1. Discretization and convergence. We begin by truncating the 
system. Consider the equation 



u t = u xx + g K (x,u,u x ), (35) 

g(x,u,u x ) for \u\ + \u x \ < K 

inf| u |+| Ua; |>x \g(x,u,u x )\ for juj + \u x \ > K 



where 

g K (x,u,u x ) := 
Consequently, 

\g K (x,u,u x )\ < \g(x,u,u x )\, 

for all x G S 1 , u,u x G K. Thanks to this, the estimates from Appendix FBI 
hold with the same constants: any complete uniformly bounded solution 



u. 



K 



(t, x) to Eqn. (J33j) satisfies 



K I + Kl + K»\ + \< I ^ c ^ W u \W> ( 36 ) 

with C independent of the truncation domain K. By choosing K appro- 
priately, solutions of Eqn. (j33|) are also solutions of Eqn. (|34|). Indeed, if 
u K (t,x) is a solution of Eqn. ()35|1 with \u K (t,x)\ < C\, then by Eqn. (jHHJ), 
\u x (t,x)\ < C2(£,Ci). If we choose K > max(Ci,C2), then solutions u K 
of Eqn. with ^^(t, x)| < Ci, are also solutions of Eqn. 
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For convenience of notation we now omit the superscript K. We discretize 
Eqn. (|35J) as follows: Let = u(t,i/d) and 

,2, ■ »'* 



Ui = d (u i+ i - 2u { + Ui-i) + g K ( -, Ui, d(u i+1 - m)j + S { (ik-i, u u u i+1 ), 

(37) 

where u[ denotes ^u(t,i/d). As before, \£f \ < |ej(d)| < C/d. The pertur- 
bations are chosen such that the given stationary solutions of Eqn. (J3"4"j) 
are also discretized solutions of Eqn. (|3*7|. 

Let be a sequence of solutions to Eqn. (|37jl with "uf(t) < Ci for 

all z and d. We will show that one can pass to the limit as d — > oo and 
obtain a complete solution to Eqn. (|34|) . The following lemma is proved in 
a manner analogous to that of Lemma ESI of 33 

Lemma 42. 

( j2^\ Au i\ 2 dt < 

J i 

where J denotes the time interval [T,T + 1], and C is independent of K. 

PROOF. If we multiply Eqn. (|37|) by Ui and then sum over i — 0, .., d and 
integrate over t G [T, T + 1] we obtain the desired estimate as in Appendix 
IBl This uses the growth of g in u x given by Hypothesis (f2). □ 



Fix K > max (C*i, C 2 ), with C\ and C 2 as above, and let /, = gx + £ 
Then 



/£i|/.l^<a 



Write each solution Ui(t) as a sum of terms Ui = u\ + uf, where 
^«f-AV = 0, ^(T) = u 4 (T), 

- AX = f h <(T) = 0. 
Then, for the homogeneous solutions uf, one estimates 

/, \ ? 1 AM|2 " 7 5 |M ' (T)|2 " C ' J ' = [T + 5 ' T + 1] - 

J i i 

This leads to the following estimate 
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For the particular solution v%, we have 

\ E I/-I 2 = i E KO'I 2 - s EW^ + \ E I^VI'- 

iii i 

For the middle term on the right hand side we have the identity — | ^(^D'^ 2 ^ = 
Jt Si ^l^ M f | 2 - Upon integration over J = [T, T + 1] we obtain 

/iE>fi 2 * = 5E|Mi? +1 f = iEi^ + Di 2 >o. 

i i i 

Combining these, we obtain 

/J E i(«n'i 2 *+ /Jei a vi 2 ^< j}j:m 2 dt<c. 

Combining the latter with the similar estimate for u\ gives the following 
estimate for the sum Ui = u\ + u\: 

/jEiwi 2 *+/jEl^.r*<c. 

Introduce the spline interpolation 

SP(wj) = £ZA 2 u i+ i(z - i/d) 3 - A 2 u i+1 (x - i/d) 2 
+ Aui(x — i/d) + Ui. 

Now set U d = SP(tij), and £/ d = PL(iij) = A-u^rr — i/d)+Ui. Then, 
\U d -U d \ 2 dxdt < ^j^O, asd^oo, 

~ C 
|17* - U d \ 2 dxdt < asrf^oo, 




J' ./s 1 




J' 



2 |2 

Mi. , 



d> 

From the latter two inequalities we derive that 

U d e H 1 ' 2 ^'- L 2 (S 1 )) n L 2 (J'; H^iS 1 )) C C(J'; H 1 ' 2 ^ 1 )), 
which implies that £\ ±|Aw;(t)| 2 < C Vt G R. Moreover, 

t/ d ,£/ d -> u, inL^J';// 1 - 2 ^ 1 )), 
L?,C/ t d - u, inL 2 (J';L 2 (^)). 
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From these embeddings one easily deduces that 

g K (x,U d ,U*) ^g K (x,u,u x ), in L 2 (J'; L 2 (S' 1 )). 



Choose smooth test functions of the form 0(t, x) = Y2k=i a k(t)wk(x), where 
{wk} is an orthonormal basis for iJ 1 ' 2 (S' 1 ). Set 4>iit) = <p(t,i/d), and 
$ d = PL(0i), then 

J ^ \ (dKji/ d, u h Auj) + £f] fadt — > J J g K (x,u,u x )(f) dxdt. 

Because of the PL approximation the following integrals become sums over 
the anchor points: 

u^ x dx = j2^ AuiA ^ = -J2^ a2u ^ 

i i 
° i i 

The final terms of the last two equations admit the following bounds: 



± E I(„; +1 -„<)A* 



< 



3d 
C 



< 7 ^o. 



lw'l 2 dx 



C 

Weak convergence implies that as d — > oo, 




J JS 1 




J JS 1 



/E^IAAl'd. 




Ufa dxdt 



j Js 1 




ucj) + | okcefc, 
w t (ia; cit, 
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where u(t,x) is the weak limit of U d (t,x). Hence, u is a weak solution 
to Eqn. for all smooth test function defined above. These func- 
tions form a dense subset in if 1,2 (J' x S 1 ), and therefore, since Ui satisfies 
Eqn. (jS3), 



Standard regularity theory arguments then yield strong solutions to Eqn. (J3"5j) 
The using the L°°-bounds on u we also conclude that u is a weak solution 
to Eqn. (|34j) . Using standard regularity techniques one can show that the 
convergence is in H l,2 (J' x S 1 ). This completes the proof of the following 
theorem: 

Theorem 43. For any sequence of bounded solutions {uf(t)} of Eqn. (J3"7j) 
with \uf{t)\ < C , for all t and i, PL(wf) converges, in H 1,2 (J x S 1 ), to a 
(strong) solution u of Eqn. (|35|). Moreover if K is chosen large enough 
then u is a (strong) solution of Eqn. (J54*j) . 

8.2. Proof of Theorem Let {u rel v} be a braid class that does 
not permit stationary solutions for Eqn. ()34j) . For d large enough the same 
holds for Eqn. (|3"Tj) : otherwise, the results in ^3 would yield stationary 
solutions of Eqn. (JMJ), a contradiction. If {u rel v} is bounded and 
proper with H(u rel v) ^ 0, then for each d large enough there exists a 
periodic solution u d with strands uf ,d (t) via j!2l Thm. 2]. By Theorem l4~3l 
this sequence yields a solution u(t,x) of Eqn. (JMJ). 

It remains to be shown that u(t,x) is periodic in t. This follows from 
the celebrated Poincare-Bendixson Theorem for scalar parabolic equations 
due to Fiedler and Mallet-Paret 11 P.. which states that a bounded solution 
u(t, x) has forward limit set either a stationary point or a time-periodic 
orbit. By assumption {u REL v} contains no stationary points which 
leaves the second option; a periodic solution. This also proves then that 
{u rel v} contains a periodic solution of the desired braid class. □ 

We remark that the proof above is for braid classes {u rel v} for which u 
has a single component. For u with multiple components, a nonvanishing 
index implies that each component of u is either stationary or periodic; 
however, unless the periods are rationally related, the entire braid class 
will be merely quasi-periodic as opposed to periodic. 
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Boundary conditions. We have employed periodic boundary conditions for 
convenience and as a means to allow for time-periodic orbits. Nothing 
prevents us from using other boundary conditions, although the resulting 
dynamics is often gradient-like. Neumann, Dirichlet, or (nonlinear) combi- 
nations of the two are imposed by choosing closed subsets B C {(0, u, u x )} 
and Bi C {(l,u, u x )} and requiring the braid endpoints to remain in these 
subspaces. As the topology of the configuration spaces of braids may 
change, so may the resulting invariants. Since the comparison principle 
still holds, our topological methods remain valid, though the invariants 
themselves may change. 

Coercivity and unbounded classes. Theorem deals with dissipative sys- 
tems. The opposite of dissipative is the coercive condition: 

uf(x, u, 0, 0) — > oo, as \u\ — > oo, 

for all x G S 1 . For either of these cases, the restriction to bounded braid 
classes may be relaxed. For dissipative systems, any braid class becomes 
bounded by adding two unlinked strands as per Appendix O In order 
to deal with coercive systems one needs to include the behavior of the 
system at infinity. We propose that a compactification of the unbounded 
braid classes yields an index with the same properties as that for bounded 
classes. 

Improper braids. A braid class is improper if components of the braid can 
be collapsed. Our results on t-periodic solutions in ^Hl dealt with improper 
braids in an ad hoc manner by 'blowing up' the collapsible strands via 
adding additional strands to the skeleton. 

A different approach would be to blow up the vector field in the tradi- 
tional manner via homogeneous coordinates, working in the setting of 
finite-dimensional PRRs. Stabilization then allows one to define the in- 
variant in the continuous limit. This type of blow-up procedure is very 
general and should be applicable to a wide variety of systems. 

Periodic skeleta. The forcing theory we have developed uses stationary 
solutions for the skeleton. We believe that all of the results hold for skeleta 
composed of time-periodic orbits. 
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p-Laplacians and degenerate parabolic equations. The fully nonlinear para- 
bolic equations studied in this paper are restricted by the 'uniform parabol- 
icity' hypothesis given by (fl). We choose to restrict ourselves to uniform 
parabolic equations in order to keep technicalities to a minimum. However, 
the theory should also apply to degenerate parabolic equations of various 
kinds. One weakening of Hypothesis (fl) would read 

< d w f(x,u,v,w), for all w^O, and (x, u, v) E S 1 x R 2 . 

Good examples of degenerate equations are the 1-dimensional porous medium 
equation u t = (u p u x ) x + g(x,u,u x ), or the p-Laplacian equation u t = 
(\u x \ p ~ 1 Ux)x + g(x, u, u x ). Solutions of these equations have less regularity 
than Eqn. (JTJ), which complicates the approach used in £0 In that case, 
one can use the weak solution approach as carried out in the periodic case. 
The key point is to find the appropriate estimates in u x . 

Scalar hyperbolic conservation laws. Conservation laws of the form 

u t = f(x,u,u x ), (38) 

where / is monotonically increasing in u x , discretize to one-sided parabolic 
systems of the form u[ = 1Zi(ui, Ui+i), cf. [15 . Our theory remains valid 
for discretized systems of this form; if we establish the appropriate a priori 
L°°-estimates a braid-forcing theory for Eqn. (|38|) can be derived. 
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The following estimates, though necessary, are antithetical to our philoso- 
phy: the entire forcing theory for Eqn. f^Q) is topological in nature. 



Appendix A. Estimates: stationary 

A stationary solution of Eqn. ((TJ) is some u 6 C 2 (R/£Z) satisfying f(x, u, u x , 
0. Hypotheses (fl)-(f2) permit the following regularity statement. 

Lemma 44. Let u e C 2 (R/£Z) be a stationary solution of Eqn. (JTJ) with f 
satisfying (fl)-(f2). There exists a constant C = C(£, \\u\\ Loo ) depending 
only on the sup-norm of u, such that 

\u x \ + \u xx \ + \u xxx \ < C. (39) 

Proof. Using (fl) we obtain the following estimate for /; 

a-(u xx )u xx + f(x,u,u x ,0) < f(x,u,u x ,u xx ) < a + {u xx )u xx + f(x,u,u x ,0), 

(40) 
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where a_ and a + are defined in |J3 Multiply Eqn. (pfUj) by u. Integrating 
over S 1 := IR/£Z, using Hypothesis (f2) and the fact that — < A" 1 yields 

u 2 x dx < / A _1 |-u| • \f(x, u, u Xl 0)\dx 
Js 

< C\\u\\ Loa / \f(x, u, u x , 0)\dx 
Js 

< C[l + [ \u x Pdx 



s 



Since 7 < 2, it follows that J s \u x \ 2 dx < C. Next we deduce from Eqn. (|3Uj) 
that \u xx \ < A _1 |/(:r, u, u x , 0)|. Again by using Hypothesis (f2) we obtain 

1 f 1 

u xx \-' = C / \f(x,u,u x ,0)\-ydx 



< 



C + C\u x \" 



dx 



< C 1+ \u r \ 2 dx < C 



The latter implies that Hull 9 2 < C . From the Sobolev embeddines for 
H /2 '^(S') we derive 

1 1 '"1 1 01,0(5) < C||«ll w 2 a < C, 

with < a < 1 — i < 1. In particular Hu^Hi, 00 < C- Again by using the 
pointwise bound \u xx \ < A _1 |/(x, u, u x , 0)| we obtain 

SUplUaa-l < C\\ f(x, U, U x , 0) || L <» 
a; 

^ C "I - C I I U-j; | | ^OO ^ Cj 

which implies that \\u xx \\ < C. By differentiating the equation and using 
the fact that / G C 1 to estimate we obtain 

dxf + d u f -u x + d u J ■ u xx + d Ux J ■ u xxx = 0. 

For u xxx this yields 

< ^{l^/l + ia^/iKi + l^/lw} < c, 



u., 



since all derivatives of / can be bounded in terms of ||u|| ioo . This completes 
the proof. □ 
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Appendix B. Estimates: non-stationary 

We repeat the regularity arguments for non- stationary solutions to Eqn. (J3J. 
As the estimates are similar in spirit as those of Appendix A, we omit the 
more unseemly steps. 

Lemma 45. Let u e C 1 (M; C 2 (M/£Z)) be a complete bounded solution with 
g satisfying (fl)-(f2). There exists a constant C = C(£, \\u\\ Loa ) depending 
only on the sup-norm ofu(t,x), such that 

\u x \ + \u xx \ + \u xxx \ + \u t \ < C. (41) 

Proof. As before, let S 1 : = R/£Z. Denote by J the time interval 
J := [T, T + 1]. Multiplying Eqn. (JSJ) by u and integrating by parts yields 

u t udxdt = — I I u 2 x dxdt+ I I g(x,u,u x )udx dt. 
jjs 1 JjJs 1 JjJs 1 

Using hypothesis (f2) we derive 

If 1 ■ 1 

u 2 dxdt<—- I u 2 dx + C + C I I \u x P dxdt. 






>j Js 1 2 j s i T j j j s i 

Hence, since 7 < 2, Jj f s i \u x \ 2 dxdt < C. 

We proceed with the more technical estimates. Given the solution u(t,x), 

2 2 , 

u t - u xx = g(x, u(t, x),u x (t, x)) EL~<(J;L~<(S )), 
since \g\ 2 ^' < C + CIm^ 2 . As such, L p regularity theory implies (see, e.g., 

INI 1 a < C(5)\\f\\ 2 2 , 

IIWsRdII 2 2 < C(5)||/|| 2 

where J' := [T + 5, T] C J for some < S <C 1. In particular, 

u E H 2 '^{S 1 )) n L°°(M; L 00 ^ 1 )). 

Bootstrapping proceeds in a standard fashion using a parabolic version 
of the Gagliardo-Nirenberg interpolation inequalities. Given any function 
u E LP(J', H^iS 1 )) n L°°(J', L^iS 1 )), then 

1 1 

\\ U \\l 2 p{J';H 1 ^p(S 1 )) — ^\\ U \\lp(J',H 2 'P) ' ll W lli°°(J'^«') > 

4^4 44 

Therefore, we have u G L~(J f ,H '"(S^)) and, hence, # G Lr 1 (J f ; L-f 7 (S )) . 
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We repeat the procedure k times, each time restricting the time domain 
[T + kS,T + 1]. Choose k > sufficiently large so that {2/j) k > 2 and 
choose 6 sufficiently small so that [T+k5, T+l] contains J" := [T+|, T+l]. 
Then we have 

u G H X > 2 {J"- L^S 1 )) n L 2 (J"; H^iS 1 )). 

By Sobolev embedding, we get u G C(J";H 1 ' 2 (S 1 )). Repeating the entire 
procedure yields u G C a { J"; C 1,Q (5' 1 )). This bound is now independent 
of T, and one translates to obtain u G C a (IR; C l,a (S 1 )). The additional 
smoothness now follows directly from the fact that u solves Eqn. 

The C 3 -estimate is obtained as in the stationary case by differentiating the 
equation and using the C^-estimates obtained above. □ 

Appendix C. Discrete enclosure 

Using a discrete version of enclosure between sub/super solutions and a 
nontrivial braid diagram, we obtain the following existence result. 

Lemma 46. Let f satisfy Hypotheses (fl)-(f2) and let v be a non-trivially 
braided stationary braid for Eqn. (jlj. Assume that there exists a u* such 
that v a (x) < u* for all a and f(x, it*, 0,0) < 0. Then, there exists a 
1-periodic solution u with 

maxi; a (x) < u(x) < u* , 

a 

for all x G S . 

It is clear that the result holds for case of a u* such that u* < v a (x) for all 
a and f(x, u*, 0, 0) > 0. In that case one finds a solution u satisfying 

u* < u(x) < minv a (x) 

a 

for all x G S 1 . 

Proof. As in $3 we discretize Eqn. (0) in x. For u* this implies 
that f(i/d, u*, 0, 0) < 0. As for the braid v we use Lemma l34l to find 
£f and the recurrence relation TZf (ui-i, U{, «i+i) := f(i/d,Ui,Aui,A 2 Ui) + 
£f(ui-i,Ui,Ui + i). By construction the discretized skeleton Disc^v is sta- 
tionary for 1Z. 

Define the region 

D = {{ui}f =0 | maxff <Ui< u*, u = u d }. 
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If the discretization is chosen fine enough then the discretized braid is non- 
trivial. As a consequence Ui cannot collapse onto Disc^v and if Ui = vf for 
some i and some a, then lZf(ui-i,Ui,Ui + i) > 0. By the definition of u* it 
follows that if Ui = u* for some i, then TZf(ui-i, U{, u i+ i) < f(i/d, Ui, 0, 0) < 
(parabolicity). The region D is therefore an attracting isolating (com- 
pact) set for Eqn. (JT5j) . Thus for each large enough d we find a discrete 
solution {uf}f =0 . Since {uf}f =0 is a priori bounded we derive from the lim- 
iting procedure in £0that this yields a stationary solution u(x) for Eqn. ((TJ), 
satisfying the desired inequality. □ 
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